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1 Introduction 


The monograph “Principles and Techniques of Shock Data Analysis” written by Kelly and Richman 
in 1969 has become a seminal reference on the shock response spectrum (SRS) [1]. Because of its clear 
physical descriptions and mathematical presentation of the SRS, it has been cited in multiple handbooks on 
the subject [2,3] and research articles [4-10]. Because of continued interest, two additional versions of the 
monograph have been published: a second edition by Scavuzzo and Pusey in 1996 [11] and a reprint of the 
original edition in 2008 [12]. The main purpose of this note is to correct several typographical errors in the 
manuscript’s presentation of a recursive algorithm for SRS calculations. These errors are consistent across all 
three editions of the monograph. The secondary purpose of this note is to present a Matlab implementation 
of the corrected algorithm. 


2 Continuous-time solution 


The shock response spectrum considers the response of a single clegree-of-freedom damped mechanical 
oscillator to a base excitation. The motion of the oscillator is y(t), and the motion of the base is x(t). This 
system obeys the following equation of motion. 

V + 20n (y - X) + ul (y - x) = 0 (1) 


Here, uj n is the natural frequency of the oscillator, and Q is the damping ratio. Accelerometer data for the 
base excitation, however, provides knowledge of x not x or x. Therefore, it is convenient to consider the 
relative motion is £(t) = y — x and the relative equation of motion. 

i + 2 (u) n £ + = — x ( 2 ) 


For arbitrary initial conditions, £o and £o> Eq. (2) has the following solution. 
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Here, Ud = o> n ^/l — C 2 - Differentiation gives the solution for the relative velocity. 
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Finally, the absolute acceleration of the oscillator can be reconstructed using the relative position and velocity. 


y(t) = C + x = -2C w n £(t) - w^C(t) 


( 5 ) 


Equations (3) and (4) correspond Eq. (4.21) and (4.27), respectively, in Kelly and Richman; however, they 
reflect several errors in the original source. 


3 Discrete-data approximation 


In practice, evaluating Eqs. (3) and (4) requires several considerations. First, the input accelerations 
are often only available at discrete time instants, and the output accelerations may only be desired at the 
same instants. The data points are here considered to be separated in time by a regular interval At. Next, 
to avoid evaluating the integrals in Equations (3) and (4) over r = 0 — > t for each instant in time, a 
recursive algorithm is desired. The solutions for Cfc = C (tfc) and Cfc = C(^fc) are treated as initial conditions 
in calculating = £(tk+ 1 ) and = £(ffc+i), and the integrals only need to be evaluated over r = 0 — > At. 
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In order to evaluate the integral terms in Eqs. (6) and (7) a continuous approximation of the discrete 
data needs to be formed. Kelly and Richman considered a parabolic approximation using the (fc— 1), k, and 
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(k + 1) data points. 
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Given a discrete time-history of accelerometer data, the resulting oscillator acceleration at any instant in 
time can now be calculated by substituting Eq. (8) into Eqs. (6) and (7), and using the resulting values in 
Eq. (5). Evaluating Eqs. (6) and (7) requires the calculation of the following Duhamel integrals. 
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Using Eqs. (9-14), the solutions for £,k+i and £,k+i can now be rewritten as shown. 
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Here, the B coefficients are defined using Eqs. (6) and (7) and the definitions of I\ through I§. 
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Finally, the expressions for these coefficients can be simplified using Eqs. (9-14). 
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(29) 
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Equations (27-36) correspond with Eqs. (6.58-67) in Kelly and Richman and Eqs. (7.38-47) in Scavuzzo and 
Pusey; however, they reflect negative-sign errors in S 3 , S 4 , and S 10 in the original source. The corrected SRS 
algorithm for a sequence of input accelerations and a desired range of natural frequencies can be evaluated 
using Eqs. (5), (15), (16), and (27-36). 


4 Comparison and Results 

The corrections described in the preceding paragraphs were verified by comparing the corrected algorithm 
and the original algorithm to an independent SRS code. The independent code used a piecewise-linear 
approximation for the base acceleration, as described by Paz [13]. The various algorithms were applied 
to accelerometer data from the ignition environment of live-fire testing of the Space Shuttle Reusable Solid 
Rocket Motor (RSRM). Specifically, data was evaluated from the radial channel at station 1479.5 on Technical 
Evaluation Motor 13. The data was sampled at 10,000 Hz. The acceleration time history is shown in Fig. (1). 
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Fig. 1: Time history of sample input acceleration. 

The SRS of this acceleration data is shown in Fig. (2) as calculated using three different algorithms. 
The corrected algorithm of Eqs. (5), (15), (16), and (27-36) are compared with the uncorrected equa- 
tions from Kelly and Richman as well as the independent code. For each algorithm, a damping ratio of 
( = 0.05 was used, and the peak response was calculated for a range of natural frequencies at one-third 
octaves up to the Nyquist frequency. The corrected algorithm and the independent code show strong 
agreement with each other; however, the uncorrected algorithm displays large differences in the high- 
frequency regime. The MATLAB implementation of the corrected algorithm is available for download at 
http:/ /www. eng. auburn.edu/users/sinclaj/SRS/. 
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Fig. 2: SRS of input acceleration for £ = 0.05 using three different algorithms. 

References 

[1] Kelly, R. D. and Richman, G., Principles and Techniques of Shock Data Analysis , The Shock and 
Vibration Information Center, Washington, D.C., 1969. 

[2] Himelblau, H., Piersol, A. G., Wise, J. H., and Grundvig, M. R., Handbook for Dynamic Data Acquisition 
and Analysis, chap. 5, Institute of Environmental Sciences and Technology, Rolling Meadows, Illinois, 
1994. 

[3] Rubin, S., “Concepts in Shock Data Analysis,” Shock and Vibration Handbook , edited by C. M. Harris, 
chap. 23, McGraw-Hill, New York, 4th ed., 1996. 

[4] Pusey, H. C., “An Historical View of Dynamic Testing,” The Journal of Environmental Sciences, Vol. 20, 
No. 5, 1977, pp. 9-14. 

[5] Houghton, J. R., “Shock Spectrum Ratios Applied to the Comparison of Pulse Signatures,” Journal of 
Mechanical Design, Vol. 102, No. 1, 1980, pp. 64-76. 

[6] Smallwood, D. O., “An Improved Recursive Formula for Calculating Shock Response Spectra,” Shock 
and Vibration Bulletin: Proceedings of the 51st Symposium on Shock and Vibration, San Diego, Cali- 




fornia, 1980. 


[7] Peleg, K. and Hinga, S., “Combining Fourier and Shock Spectra for Improved Test Specification,” 
Journal of Sound and Vibration, Vol. 110, No. 2, 1986, pp. 289-307. 

[8] Smith, S. and Melander, R., “Why Shock Measurements Performed at Different Facilities Don’t Agree,” 
Proceedings of the 66th Shock and Vibration Symposium, Biloxi, Mississippi, 1995. 

[9] Sek, M., Rouillard, V., Tarash, H., and Crawford, S., “Enhancement of Cushioning Performance with 
Paperboard Crumple Inserts,” Packaging Technology and Science, Vol. 18, 2005, pp. 273-278. 

[10] Shin, K. and Brennan, M. J., “Two simple methods to suppress the residual vibrations of a translating 
or rotating flexible cantilever beam,” Journal of Sound and Vibration, Vol. 312, 2008, pp. 140-150. 

[11] Scavuzzo, R. J. and Pusey, H., Principles and Techniques of Shock Data Analysis, The Shock and 
Vibration Information Analysis Center, Arvonia, Virginia, 2nd ed., 1996. 

[12] Kelly, R. D. and Richman, G., Shock Data Analysis - Engineering Principles and Techniques, Wexford 
Press, Eldersburg, Maryland, 2008. 

[13] Paz, M., Structural Dynamics: Theory and Computation, Van Nostranc! Reinhold, New York, 2nd ed., 
1985, Sect. 4.4. 


9 



